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ABSTRACT 


Transport imaging is a SEM-based technique used to 
directly image the motion and recombination of charge in 
luminescent semiconductors, allowing for the extraction of 
transport parameters critical to device operation. In this 
thesis, transport imaging for ID structures was initiated 
with work on sample preparation, modeling and initial 
characterization. One dimensional structures are being 
integrated into forefront electronics due to their inherent 
advantages in size, packing density and power consumption. 

In this work the one dimensional equation for steady state 
minority carrier recombination distribution solved for the 
Gaussian source is derived and results from numerical 
simulations are presented. The diameter of the SEM beam is 
determined experimentally allowing for accurate simulation 
parameters. Intensity and drift measurements on four 
batches of top-down wire structure samples, fabricated on a 
AlGaAs/GaAs/AlGaAs double heterostructure using a EIB, are 
presented. Significant decreases in luminescence in EIB 
exposed regions are reported. Spatial luminescence from 
single bottom-up GaN and ZnO nanowires deposited by metal 
initiated metal-organic CVD on Au and Si02 substrates is 
imaged. CL spectra for GaN and ZnO, with peak intensities 
at 3.27 and 3.29 eV, are characterized. Einally, several 
suggestions for further research are offered including 
transport imaging on contacted bottom-up nanowires and a 
potential application of transport imaging to EIB damage 
characterization. 
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I. 


INTRODUCTION 


A. BACKGROUND 

The transport of free charge in semiconductor 
structures is fundamental to the operation of many modern 
optoelectronic and semiconductor devices. Photodiodes 
including solar cells, photomultipliers used in night 
vision systems and spontaneous and stimulated emission in 
Light Emitting Diodes (LEDs) and lasers are just a few of 
devices and physical processes in which free charge 
transport parameters such as carrier mobility are of 
critical importance. 

Since the early 1970's much research in optoelectronic 
and semiconductor devices has been focused on decreasing 
device sizes due to both practical advantages, including 
increased packing densities and reductions in device 
production costs and power consumption, as well as physical 
advantages such as dimensional differences in the density 
of states (DOS)[1]. Quantum wells, fundamentally two 
dimensional (2D) structures, have found myriad applications 
in the fabrication of laser diodes, infrared (IR) quantum 
cascade lasers and High Electron Mobility Transistors 
(HEMTs) to name just a few [2] . Research in the field of 
one dimensional (ID) structures, referred to as nanowires, 
has already yielded many advances including nanowire field 
effect transistors (FETs) , crossed nanowire p-n junctions, 
InP nanowire LEDs and even ultraviolet (UV) lasers made of 
single GaN nanowires much like the one pictured in Figure 1 
[3,4,5]. Accurate measurement of transport parameters is an 
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important part of nanowire applications for all the above 
mentioned nanowire devices. 



Figure 1. 7500x SEM Micrograph of GaN nanowire on Si02 

Substrate 


Traditional measurement of transport parameters in 
materials is done electrically, as in diffusion length 
measurements performed through electron beam induced 
current (EBIC) experiments. Electrical measurements are 
often affected by contact related phenomena that require 
potentially expensive and damaging sample processing. In 
view of this, a contact free technique called 'transport 
imaging' was developed in our laboratory in 2004 [6]. 

Transport imaging employs an optical microscope (OM) , 
internal to a scanning electron microscope (SEM), in order 
to directly image the spatial recombination of charge 
generated by electron irradiation at a point. Planar 
contacts can be used to create local fields so that both 
diffusive and drift behavior can be directly observed. 
This has been successfully performed in two dimensional 

2 




heterostructures [6,7] and in bulk materials. One goal of 
this thesis is to advance the technique by investigating 
its initial applications to ID structures. 

The SEM-based approach also provides the benefit of 
high resolution for excitation and an ability to easily 
excite wide bandgap materials. The transport imaging method 
makes it possible to extract key materials parameters, such 
as diffusion lengths and minority carrier mobility, from a 
single charge coupled device (CCD) picture. 

B. MILITARY RELEVANCE 

Solar cells and integrated circuits (ICs) on 
intelligence and communication satellites and laser target 
designators and night vision goggles (NVGs) in the 
battlefield are just a few of the roles optoelectronic and 
semiconductor devices play in the modern military. Their 
importance cannot be overstated, and neither can the 
importance of the development of physical understanding and 
processes to advance these technologies. 

One of the next steps in many of these technologies is 
the fabrication of nanowire based devices. A single 
nanowire, for example, can function as a stand alone cavity 
and gain medium, acting as a very narrow linewidth emission 
laser which could have many applications, including optical 
computing and microanalysis [4]. Reduced lasing threshold 
currents (see Figure 2) and diminished temperature 
dependence of the threshold current observed in nanowire 
lasers mean more efficient and compact uncooled laser 
devices [5]. Higher modulation factors in electro-optical 
modulators and new nanowire FETs may allow for faster, more 
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compact, efficient and inexpensive optical integrated 
circuits (OICs) and ICs [2] . As a result, DARPA and other 
agencies that support research and development for military 
optoelectronics are investing significant resources in 
nanotechnology. 



Figure 2. Evolution of Threshold Current for Semiconductor 

Lasers (adapted from [2]) 

Figure 2 shows the impact of the double 
heterostructure on threshold currents for semiconductor 
lasers. The result of the movement from 3D bulk materials 
to 2D double heterostructures allowed for reduction in 
threshold currents, lower power consumption and a reduction 
in physical dimensions that has put double heterostructure 
semiconductor lasers in every compact disk reader and many 
other optoelectronic devices. The military applications of 
nanowire technologies in semiconductor and optoelectronic 
devices could result in the same widespread revolution in 
technology. 
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C. THESIS OVERVIEW 

In this work the transport imaging technique will be 
expanded to examine spatial luminescence from ID structures 
in order to develop a contact free method of measuring 
carrier transport parameters. Transport is a key 
application of quantum wires. The possibility of direct 
measurement allows for more accurate transport parameter 
measurements and visualization of spatial variations along 
the length of the nanostructure. Effective application of 
the transport imaging method should allow the determination 
of values for minority carrier diffusion length and 
potentially, minority carrier mobility. These are important 
for devices, such as lasers, that depend on minority 
carrier injection. Transport characterization in quantum 
wires is still in a relatively early stage and the 

transport imaging technique would provide a unique way to 
extract important information. 

Chapter II discusses the theoretical and experimental 
determination of the size of the SEM beam as well as the 

dynamics of electron beam-material interaction and carrier 
generation, critical to the SEM-based study of nanoscale 
structures. The one dimensional equation for steady state 
minority carrier recombination distribution solved for the 

Gaussian source is derived and results from numerical 

simulations based upon the previously developed equation 
are presented. 

In Chapter III the experimental apparatus as well as 
the fabrication of multiple batches of wire structures on 
Beryllium (Be) modulation-doped p-type AlGaAs/GaAs/AlGaAs 
heterostructure using focused ion beam (FIB) lithography 
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are described. GaN and ZnO quantum wires grown by chemical 
vapor deposition (CVD) are also described in detail. 

Chapter IV present the results of the luminous 
intensity and drift experiments carried out on the four 
separate batches of top-down wire structures and the 
cathodoluminescence (CL) spectra and results of different 
imaging techniques on the bottom-up nanowires. The final 
chapter provides an analysis of the previously obtained 
results and presents conclusions for both the top-down and 
bottom-up wire structures as well as ideas for future 
research. 
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II. TRANSPORT IMAGING THEORY AND MODELLING 


A. INTRODUCTION 

Transport imaging utilizes the combination of a SEM 
for high resolution charge generation and a Silicon CCD in 
conjunction with an OM for the recording of spatially 
resolved luminescence. Upon exposure to the electron beam, 
a steady state spatial distribution of luminescence 
associated with charge recombination is created which can 
then be recorded by the CCD. The CCD images can be analyzed 
in order to provide quantitative measurements of local 
minority carrier diffusion and drift lengths. 

In order to extract minority carrier transport 
parameters from CCD images, the electron beam- material 
interaction and the subsequent generation, drift, diffusion 
and radiative recombination of carriers in the material 
must be modeled. The SEM beam's effects on the sample can 
be modeled as an extended generation region within the 
optically active GaAs layer of the sample. By employing 
certain simplifying assumptions, a model can be developed 
from the continuity equation which directly relates the 
spatial luminescence associated with minority carrier 
recombination to the minority carrier's diffusion and 
drift. 

B. MODEL PARAMETERS AND ASSUMPTIONS 

1. Gaussian Generation Region and Beam Size 

It has been established experimentally that the SEM 
beam's intensity follows a Gaussian distribution to first 
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order. Theoretical calculations of the beam size assume 
that the beam diameter, defined as the Full Width at Half 
Maximum (FWHM) of the spot distribution, is the quadrature 
sum of the Gaussian probe electron beam and the various 
aberrations and has a beam diameter of ~50 nm at 30 keV and 
lxl0~^° Amperes (A) [8] . 

In order to model the generation region, the goal was 
to experimentally determine the beam diameter of our SEM, 
because of potential deviation from theoretical or 
published values during operating conditions. A Geller 
Microanalytical Laboratories Magnification Reference 
Standard (MRS-3), pictured in Figure 3 below, was used to 
provide a high contrast, high resolution pattern. 



Figure 3. 20x SEM Micrograph of the MRS-3 Sample 

Using the MRS-3, the rate of change of intensity along 
sharp edges of the control sample was correlated to the 
diameter of the SEM beam experimentally. The SEM beam 
profile over a discontinuity was measured under a constant 
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accelerating voltage of 30 keV at different probe currents, 
as shown in Figure 4. 



Figure 4. Probe Current Dependence of 30 keV SEM Beam Profile 

Over Sample Discontinuity 

The beam intensity profiles obtained experimentally 
were analyzed and an empirical value for the SEM's beam 
width was calculated as a function of probe current. Values 
for 2a, where <7 is the standard deviation of the Gaussian 
distribution, were obtained from the positions of the 16% 
and 84% thresholds for beam intensity and converted into 

beamwidths following FWHM = 2iy2ln(2)cr « 2.35cr [ 9 ] . 


Probe Current 

16%-84% (2a) 

FWHM 

(A) 

(nm) 

(nm) 

1x10“^^ 

34 

40 

3x10“^^ 

39 

46 

1x10“^° 

48 

56 

3x10“^° 

107 

126 

1x10'® 

139 

163 


Table 1. 30 keV SEM Beamwidth Dependence on Probe Current 
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The expected behavior as a function of probe current 
was observed as increasing the probe current increased the 
beam diameter and decreased the spatial resolution of the 
SEM. The SEM beam was measured to be a well focused beam 
(FWHM ~50 nm) at low probe currents. A similar range of 
experiments could be carried out varying the accelerating 
voltage for a given probe current. 

Modeling of the beam size's relation to the shape of 
the luminescent distribution shows a relative degree of 
insensitivity, especially when the beam FWHM << Ld, the 
diffusion length. This is even the case in the ID diffusion 
scheme, which characteristically shows increased 
sensitivity to generation region variations. Figure 5 shows 
a comparison between ID generation regions of 28 nm (n=10) 
and 166 nm (n=60), for an Ldof 3.6 ym. 



Figure 5. ID Beamwidth Comparison, L^=3.6 pm 
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The Gaussian distribution function f (x) =—j= e " , where 

yjTT 

2 1 

n —- can then be used to model the generation region of 

the electron beam from the SEM, where the value for n~30 
corresponds to a 30 keV, 1x10”^° A electron beam of width 56 
nm. To first order, in the heterostructure material of 
interest, the generation region created by the incident SEM 
beam is approximated by the incident beam diameter on the 
surface, since the width of the layer is only ~1000 A. 

2. The Low Injection Limit 

In order for the spatial luminescence information 
recorded by the CCD to be a direct indication of the 
minority carrier distribution without resorting to more 
complex models, the minority carriers must be recombining 
into a much larger and effectively constant number of 
majority carriers. This is known as the low injection 
limit. 

The number of holes and electrons per unit volume, p 
and n, can be calculated in the optically active GaAs 
region of the heterostructure to be ~10^®cm”^ and ~10®cm”^ 
respectively, at room temperature. Ap and An would then be 
the number of holes and electrons per unit volume generated 
by the electron beam incident on the material. In order for 
the low injection limit assumption to hold for the 
modeling, Ap«.p while An » n . 

The maximum electron-hole pair generation rate, 
neglecting electron beam energy loss due to backscattering, 
can be calculated by G = £’^/£’,., where is the incident 
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beam energy and E- is the ionization energy of an electron- 
hole pair in GaAs. For electron beam energies in the 5-35 
keV working range of the SEM, is independent of incident 
electron beam energy and has been empirically determined to 
be 2.8 [10], where E^ is the bandgap of GaAs, ~1.42 eV at 

300 K. The total minority population is then —Gr where I is 

e 

the probe current and r is the minority carrier lifetime in 
the GaAs heterostructure, ~4 ns. 

The minority carrier concentration in the generation 
region can then be calculated by using a hemispherical 
generation volume approximation for a 30 keV, 10 nA 
electron beam, following the development of Kanaya and 
Okayama [11] . This yields a minority carrier concentration 
of ~10^®cm”^. This value is an overestimate of the minority 
carrier concentration due to surface recombination and the 
dimension of the sample (0.175 ym) relative to that of the 
beam penetration (2-5 ym) and is still over two orders of 
magnitude lower than the effective majority carrier 
population due to doping of the GaAs in the samples of 
interest. All experimental work will be done below the low 
injection limit of a 30 keV, 1 nA electron beam. 

Within the low injection limits the model is valid, 
but some effects are not physically observable. In the case 
of very short diffusion lengths in materials such as GaN 
(L~0.1 ym) that approach the optical imaging system's 
resolution of ~0.4 ym, the simulated effects cannot be 
observed. In very large diffusion length materials such as 
GalnAs solar cell heterostructures (L~130 ym) it is very 
difficult to stay within the low injection limits as the 


12 



carrier lifetimes can increase over three orders of 
magnitude, from ns to ms. In addition, recombination in 
these materials occurs over such a large area and can be 
below the required intensity per unit area to register 
above noise levels on the CCD. 

C. ID STEADY STATE MODEL 

1. The Non-homogeneous Second Order Differential 
Equation 

The top-down wire structures described in Chapter III 
vary from ID and y,z<L^) quantum wires to 2D {x,y:^L^ 

and z<L^) heterostructures in their physical dimensions. 
The development of steady state carrier recombination 
models for all the wire structures is required; however, 
the model for 2D structures has been developed previously 
[6,7] . A ID steady state equation is now developed, 
starting with the continuity equation for electrons in a p- 
type material: 

^ = G,-V+-(VU,) (1) 

at q 

where the time rate of change of electrons per unit volume, 
dn 

—, is expressed in terms of and U^, the generation and 
dt 

recombination rates per unit volume, and the divergence of 
J^, the current density vector. In the low injection limit 
An 

U ^=— where An and are the number of minority carriers 
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per unit volume and their respective lifetime. The current 
density for minority carriers can be defined as: 


Jn = + qDyn 


( 2 ) 


where is the mobility of electrons in GaAs, E is the 

applied electric field within the material and is the 

diffusion coefficient for electrons. Substituting Equation 
(1) into Equation (2) we arrive at: 


dn ^ An 1 

-+- 

dt r„ q 


V • {q/u„nE + 


(3: 


expression, can 

be related to 

the 

L, 

by E)^= —, An»n 

since and. 

by 

the 

steady state 

there is no time 

= 0 . 

We also realize 

that any gradient 

or 


dependence, so — 

dt 

external electric field will only exist in one dimension, 

d^n 


arbitrarily x, such that E = qEx and V n = . Equation (3) 

can now be written as: 


dx^ 


_ ^ n ^dn L d n 

0 = G„- + MnE — + — ^ 

T„ dx r„ dx 


(4) 


Equation (4) can be further simplified by the substitution 
S = ju^T^E, where S is the drift length, and rearranged into a 
non-homogeneous second order differential equation: 
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(5) 


d^n ^ S (dn^ 


dx" 


\dxj 


1 G T 
- n =- — 

1 } L" 


2. General Solutions for a Point Source 

The derivation of the one dimensional equation for 
steady state minority carrier recombination distribution is 

first solved for the point source, such that G^=^^S(x) 

where d{x) is the delta functional. Equation (5) then 
becomes: 


1} 



^ dn^ 


-n = -gLS(x) 


( 6 ) 


The Fourier transform can be taken of Equation (6) 
yielding: 


-l}k^N{k) - ikSN{k) - N{k) = -gL 


(7) 


Where N(k) is the Fourier transform of n. Equation (7) can 
be re-arranged to form: 


Nik) = ^ 

L 


r 


yk^ + CJi + C2 j 


( 8 ) 


where the 

constants c. 

II 

and r • 

The operation 

of 

completing 

the square 

can 

then be 

performed on 

the 

denominator 

of the expression 

in Equation 

(8) 
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k^+CJi + C^= k + ^ +C3" 

V 2 J 


^ ^S^+AU ,, , C, 

where C 3 =-r-. After a change of variables k = k -we 


arrive at: 


N(k) = - 


L{(kr+c, 


-i 

Equation (10) can then be multiplied by — and inverse 

C 3 

Fourier transformed: 


n(x) = [ 

rr J 




{k'f+C, 


Equation (11) can be recognized as a Fourier transform of 
an exponential in conjunction with a multiplicative 

q. 

- -IX 

exponential term, e ^ , yielding the one dimensional 
equation for steady state minority carrier recombination 
distribution solved for the point source: 


n(x) = ■ 


sx 


_ ^ ^ 


S^+4L^ 
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3. Solution for a Gaussian Generation Region 

For the case of the source function, p (x) = r the 

T 

solution to the point source can be generalized, as long as 
the function is a general member of a delta sequence, such 

oo 

that lim \ S^(x)f (x)dx — f (0) . The Gaussian distribution function 

n—>(X) J 
—00 

n 2 1 

f{x) = —j=e , where n =—satisfies these conditions. 

V;r 2cr 

Following the development preceding Equation (10) a 
parallel equation can be derived which takes the modified 
source function into account: 


Nik) = 



(13) 


where F(k) is the Fourier transform of f (x) . The inverse 
Fourier transform of N(k), w(x), can be expressed as a 

convolution, that is = f * g ■. 


w{x) 


F(k) 


L (V) +C, 


Applying the convolution theorem, 

GO 

f*g=f fix') 


uu 

j G{k)e 


ik(x-x') 


dk 


dx' 


(14) 


(15) 
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and taking advantage of the previously derived result in 
Equation (12) for the result of the Fourier transform of 
the second term. Equation (16) is obtained: 


w(x) = j 


1 n 


lU 




S{x-x') 

9 

e e 


^/s^+4Lh 


dx' 


J 


(16) 


After rearranging. Equation (17) is final form of the one 
dimensional equation for steady state minority carrier 
recombination distribution solved for the Gaussian source: 


w(x) 


2n 




'+4L" 


S(x-x') 

, 21 } . 




+4r 


2r 


dx' 


(17) 


D. MODELING RESULTS 

1. Introduction 

The steady state carrier recombination models 
developed in the previous section can be used to calculate 
expected spatial intensity distributions. MATHCAD was used 
to calculate the single and double integrals numerically, 
employing a Romberg trapezoidal approximation with 
constraint and convergence tolerances of 0.001. Numerically 
the integration over an infinite generation region is not 
possible, so a finite integration over three standard 
deviations of the Gaussian generation region accounting for 
99.7% of the total generation is used as an approximation. 
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2. 


Non-Normalized Diffusion Modeling 


Non-normalized diffusion modeling of top-down wire 
structures shows a substantial expected difference in 
recombination intensity between the ID quantum wires and 2D 
heterostructures, assuming relevant material parameters 
stay constant. At its peak, the expected intensity 
difference is approximately 30 times greater for a ID drift 
scheme than a 2D drift scheme, as illustrated in Figure 6. 
On the logarithmic scale the slope of the ID line is 
linear, while the function of the 2D line is a zeroth order 
modified Bessel function of the second kind whose slope is 
approximately linear far from the origin. 

0.07 

0.06 

O.OS 

C 0.04 
O 
C 

o 0 03 
3 

S 0.02 
< 

0.01 

0.00 


Figure 6, ID and 2D Intensity Comparison, Linear and Log 

Scales, L^=3,6 pm, n=30 

3. Normalized Drift Modeling 

While the non-normalized results demonstrate the clear 
difference between ID and 2D intensities, normalization 
allows for a more effective comparison of the predicted 
shapes of spatial intensity distributions. This method of 
comparison is useful because of the absolute intensity 
fluctuations due to SEM optics and natural filament 
thermionic emission variations between experiments. Figures 
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Normalized Intensity Normalized Intensity 


7 and 8 show a change of three orders of magnitude in 
applied field on ID and 2D structures and the resulting 
drift on both linear and logarithmic intensity scales: 



Figure 7. ID Drift Intensity Comparison, Linear and Log 

Scales, Ld=3.6 pm, n=30 



Figure 8. 2D Drift Intensity Comparison, Linear and Log 

Scales, Ld=3.6 pm, n=30 


4. Comparison of ID and 2D Normalized Drift Modeling 

In the modeling process an arbitrary amount of 
electron-hole pairs are created in the generation region 
and their movement and recombination are simulated. By 
using the same generation regions and the same material 
parameters in the ID and 2D models, only the dimensional 
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parameters affect the intensity and distribution of the 
luminescence. 

In the 2D case, the generated charge drifts and 
diffuses into a region (2nrdr) which increases in direct 
proportional to the distance traveled. In the ID case 
however, there is no change in the size of the region (dr) 
charge diffuses into with distance traveled. In practice, 
when the carrier diffusion lengths become much greater than 
the width of a sample heterostructure, an effectively ID 
drift and diffusion scheme is created. 

A comparison between ID and 2D normalized drift 
characteristics shows differences in the shape of luminous 
intensity that increase with applied field. With increasing 
drift (ju^r^E) lengths, proportional to increasing applied 
external fields, the effective confinement experienced by 
the minority carries increases. Due to this effect, 
transport in ID is more sensitive to variations in material 
parameters. Figure 9 shows variation between ID and 2D low 
and high applied field drift conditions on both linear and 
logarithmic scales: 



Figure 9. ID and 2D Drift Intensity Comparisons, Linear and 

Log Scales, Ltj= 3.6ijm, n=30 
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III. EXPERIMENTAL APPARATUS AND SAMPLES 


A. EXPERIMENTAL APPARATUS 
1. The SEM 

A JEOL 840A SEM with an internal optical microscope 
(OM) and mirror assembly allowed for SEM, optical and CL 
imaging. An OM objective connected to a charge coupled 
device (CCD) camera or a mirror assembly connected to a CL 
system was inserted by means of retractable arms into the 
SEM chamber for data collection. The SEM system uses a 
constant flow liquid cooled stage which operates between 
room temperature and about 5K if liquid helium is used as 
the coolant. Figure 10 shows a schematic representation of 
the SEM system. 



The Oxford Instruments CL system connected to the 
retractable mirror employs a parabolic mirror, H m 

monochromator and TE cooled GaAs PMT as a detector. An 
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Apogee silicon CCD camera with a 2184x1472 pixel array 
serves as the optical detector connected to the retractable 
OM. 

The Apogee CCD camera collects unfiltered light from 
about 300-1100 nm on a 15x10 mm CCD array. The CCD camera 
is thermoelectrically cooled to temperatures below -10°C to 
reduce noise during operation. Each pixel on the array has 
dimensions of 6.8 ym x 6.8 ym and the OM has a 
magnification of approximately 20x. This yields a CCD on 
sample resolution of about 0.4 ym/pixel. The resolution for 
incoherent light can be estimated as 

0.6U 

(IS) 

NA 

where Ay is the wavelength dependent spatial resolution, A 
is the wavelength and NA is the numerical aperture. For the 
room temperature peak emission of GaAs (870 nm) and a 
numerical aperture of 0.95 (as a maximum value). Ay is 
found to be 0.56 ym. Thus the CCD resolution is close to 
the diffraction limit of 870 nm light characteristic of 
GaAs . 



2. SEM Modes 

The SEM can operate in picture, line and spot modes. 
Line and picture modes allow the beam to raster in a 
straight line or over a rectangular area, while spot mode 
focuses the beam down, as much as possible, to a point. SEM 
operating modes and magnifications can be chosen depending 
on the specific generation region shape and size desired. 
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In spot mode, the electron beam is focused down to a 
single point whose actual size depends primarily on probe 
current and accelerating voltage as discussed in Chapter 
II. Spot mode also displays a small secondary spot due to 
internal reflection in the OM optics, as seen in Figure 11. 



Figure 11. 30 keV, 1x10 A, 1 Second Spot on GaAs 

Heterostructure (154 pm x 125 pm) 

Imaging with the SEM as well as sample measurements 
can be done in line and picture modes. Picture mode also 
serves as an ideal tool for fine focus adjustments in 
height for the optical system and can serve as a method for 
identifying important sample features relative to a spot 
when used in conjunction with spot mode operation during 
the same CCD exposure. 
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3. 


Data Extraction 


Diffusion and drift spot data in the ID limit takes on 
an exponential form, as described in Chapter II, which can 
be directly extracted from the value of the slope in the 
CCD pictures intensity profiles. Data from the line scan 
feature included in the Micro CCD software package or data 
from an imported image matrix in MATLAB can be used to 
obtain intensity profiles from CCD images. To reduce noise 
and locate and extract luminescence maxima from the TIF 
files in which CCD data is recorded a simple MATLAB program 
was used (see Appendix A). 

B. SAMPLE PROPERTIES AND FABRICATION 

1. Sample Heterostructure 

All top-down samples were fabricated from a 
AlGaAs/GaAs/AlGaAs double heterostructure grown by solid 
source molecular beam epitaxy (MBE) on a semi-insulating 
substrate. Be was used to modulation dope the AlGaAs 
regions, making the optically active GaAs layer p-type. The 
minority carriers are electrons and within the GaAs layer 
there is an effective hole concentration of ~lxlO^® cm“^. 
Figure 12 illustrates the heterostructure dimensions. 

Luminescence from the semi-insulating substrate is 
negligible compared to that of the epitaxial GaAs in CL 
measurements due its significantly lower level of optical 
activity relative to the GaAs. 
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Figure 12. The AlGaAs/GaAs/AlGaAs Heterostructure 

2. Focused Ion Beam (FIB) Lithography in Wire 
Structure Fabrication 

In order to create the 'top-down' wire structures for 
carrier diffusion and recombination studies in this thesis, 
one of the methods used was FIB lithography. FIB is a 
technique in the semiconductor industry used in 
applications such as mask repair and site specific 
Transmission Electron Microscope (TEM) sample preparation. 
Exposing the heterostructure to the EIB served as a means 
of physically and electrically isolating the wire 
structure. 

Inside the EIB system, a Ga liquid metal ion source 
(LMIS) is used to create ionized Ga ions which are then 
accelerated through a potential difference of 10-40 keV and 
focused with electrostatic lenses. When the high energy Ga 
ions strike the sample, the surface is amorphized and atoms 
on the surface are sputtered, which allows the sample 
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exposure to a FIB to be analogous to milling the exposed 
areas. 

The FIB system also contains an SEM which operates at 
an offset of 52° from the FIB. The sample is connected to a 
stage which can be rotated to align it with the electron 
beam or the ion beam. Short of using the ion beam to image, 
however, the FIB etching must be done while viewed from a 
52° offset which introduces challenges to milling exact 
structures . 

3. Chemical Vapor Deposition (CVD) in Nanowire 
Fabrication 

The 'bottom-up' nanowire structures studied in this 
thesis were deposited by metal initiated metal-organic CVD 
on Au and Si02 substrates. The process was carried out with 
different reagents in order to deposit both GaN and ZnO 

nanowires on the substrate at temperature of 800~1000 °C 

[14] . CVD is widely used in the semiconductor industry for 
thin film deposition. 

In CVD the substrate is placed on a temperature 
controlled susceptor, made from a material resistant to the 
metalorganic compounds. Nitrogen gas is then bubbled 
through a metalorganic liquid and used as a carrier gas. In 
the case of GaN nanowires, nitrogen was bubbled through 
Trimethylgallium (TMG) to create a metalorganic vapor and 
introduced into the sample chamber with ammonia gas. 

Thermally evaporated metals 2-10 nm thick on the 

substrate such as Ni, Fe and Au are then used as initiators 

for vapor-liquid-solid (VLS) nanowires. The VLS process 

begins with the dissolution of gaseous metalorganic 
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reactants into the liquid droplets of the catalyst metal. 
Nucleation and growth of singlecrystalline rods can then be 
achieved as the metalorganic vapors combine in on the seed 
crystal, allowing the growth of nanowires [14]. 

C. FIRST SAMPLE BATCH 

1. Batch Fabrication 

The fabrication of the first sample batch of top-down 
structures began with a methanol rinse of the fragment of 
the AlGaAs/GaAs/AlGaAs heterostructure to be etched by the 
FIB. In order to create an electric field within the sample 
for potential drift measurements. Platinum (Pt) pads for 
contacts were made on the sample by electron beam 
deposition. 

In order to align the wire structures with the Pt 
pads, the sample was imaged by a 500 pA ion beam for a 
number of seconds before it was etched by the FIB, with 8- 
10 minute exposures to a 500 pA Ga ion beam. The electron 
beam deposition of Pt contacts and the subsequent imaging 
and etching by the FIB was repeated to fabricate all of the 
wire structures. 

The sample was then rinsed in solution of 
Hydrocholoric Acid (HCl) for 30 seconds to remove 
accumulated Ga on its surfaces produced by the FIB etching 
process. 


2. Batch Sample Dimensions and Details 
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The first sample batch was fabricated with widths 
ranging from approximately 1 to 40 ym. High resolution SEM 
micrographs revealed variations in the wire widths of up to 
±100 nm, as illustrated in Figure 13. The irregularities in 
the attempted fabrication of uniform width structures also 
brings first sample batch fabrication methods into 
question. In Table 2 the error column is a combination of 
both measurement uncertainty and structure width variation. 



Figure 13. lOOOx SEM Micrograph of 923 nm Wire Structure 


Sample 

Number 

Wire Width 

(nm) 

Error 

(nm) 

Error 

(%) 

1 

39150 

±100 

0.26 

2 

9977 

±45 

0.45 

3 

948 

±35 

3.69 

4 

923 

±40 

4.24 


Table 2. First Batch Sample Dimensions 

3. First Batch Damage Considerations 

Measurements on and near the FIB etched wire 
structures displayed luminous intensities over 1000 times 
lower than measurements taken a few millimeters away and 
served as a clear indication of material damage associated 
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with decreased minority carrier lifetime in the vicinity of 
the wires. In order to understand and quantify the decrease 
in intensity in the vicinity of the wire structures, a 
series of CL spectra were taken. Figure 14 shows the four 
FIB etched wires, indicated by white arrows and the 
position of the CL spectra taken, indicated by the numbers 
1 to 8. The spectra data points were evenly spaced in 
intervals of 465 ym. 



Figure 14. 550x SEM Micrograph of First Batch Sample CL Points. 
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Figure 15. Spatial Variation of CL Spectra on First Batch 

In the vicinity of the wire structures CL luminosity 
decreased significantly, and was only measurable in the 
areas where the substrate was exposed by the FIB process, 
as shown in Figure 16. 



Figure 16. SEM and Corresponding CL Data for 40 pm Wire 


Without the smaller bandgap heavily doped GaAs 

material covering the semi-insulating substrate in the FIB 

etched regions it was determined that observed luminosity 

came directly from the substrate rather than the epitaxial 

GaAs layer. Figure 17 shows broadening of the spectra and 

shift in energy peaks typical of heavily doped materials, 
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when comparing the doped GaAs and the FIB etched well 
substrate regions of the 40 ym wire structure. 



Figure 17. FIB Etched Well Area and Epitaxial Layer CL Spectra 

The data collected on the first sample batch was clear 
evidence of large scale damage to the heterostructure 


material 

as 

well 

as 

large wire structure 

width variation. 

Exposure 

to 

the 

FIB 

imaging 

was enough 

to amorphize and 

etch surface 

of 

the 

material 

in the vicinity of the wire 


structures, altering its optical properties significantly, 
more than expected [15]. 

D. SECOND SAMPLE BATCH 

1. Batch Fabrication 

In order to reduce the damage to the sample that 
occurred in the fabrication of the first batch, isolating 
the wire structure region from FIB damage became a 
priority. In order to achieve this, contacting on the 
sample had to be done in a way that did not require FIB 
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imaging of the area. A 200 ym Silicon mask was affixed 
across the selected AlGaAs/GaAs/AlGaAs heterostructure 
fragment with wax to protect the area on which the wires 
would be etched from electron beam deposition of Ni/Ge/Au 
contacts that followed. Figure 18 gives a schematic 
representation of the area processing of the sample batch. 



The silicon mask was then physically removed and 
methanol and xylene chemical rinses were used to clean the 
area. Rapid thermal annealing (RTA) of the contacts was 
then carried out at 430 °C for 20 seconds and the wire 
structures made in the non-metalized region of the sample 
using 8-10 minute 500 pA FIB etching. The sample batch was 
then sonicated in solution of HCl, causing the metal 
contacts to become detached from the sample. 

The contact areas were then cleaned with low power 
oxygen plasma and the Si mask was reattached to the sample 
with wax. A second attempt at electron beam deposition of 
Ni/Ge/Au contacts also failed and after sonication in 
methanol and xylene chemical rinses the sample was cut with 

a wire saw, to electrically isolate each wire structure, as 
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show in Figure 18. Pressed In contacts were added to 
individual samples as required for application of external 
bias. 


2. Batch Sample Dimensions and Details 

The second sample batch was fabricated with wire 
structure widths which ranged from about 0.7 to 24 ym. No 
visual evidence of amorphization on the sample surface due 
to exposure to the FIB was present, but partially deposited 
contacts which failed to adhere to the sample surface are 
clearly visible in Figure 19. 



Sample 

Number 

Wire Width 
(nm) 

Error 

(nm) 

Error 

(%) 

1 

23630 

±85 

0.36 

2 

13890 

±50 

0.36 

3 

1725 

±20 

1.16 

4 

746 

±35 

4.69 


Table 3. Second Batch Sample Dimensions 
Figure 19. lOOX SEM Micrograph of Contacted Second Sample Batch 


A significant increase in the consistency of the wire 
widths was also observed, as shown in Figure 20. 
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Figure 20. 


7500x SEM Micrographs of 1.7 pm and 14 pm Structures 


3. Second Batch Damage Considerations 

In the second sample batch optical measurements near 
the FIB etched wire structures displayed luminous 
intensities over 10 times higher than measurements taken in 
the same regions of the first sample batch. While less 
damage was done to the samples, the absolute intensity in 
these regions was still almost 100 times lower than that of 
the control sample. 

The 700 nm wire demonstrated both light and dark 
regions under SEM and CL measurements and visible variation 
in its width and structure. All wire structures also 
demonstrated surface irregularities in the regions 
surrounding the FIB milling. The second batch samples were 
exposed to multiple RTA and chemical processes and the 
observed data indicate that these processes could have 
damaged or degraded the materials. Figure 21 shows the 
surface irregularities in the vicinity of the 700 nm wire 
as well as the irregularities within it. 
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Figure 21. 550x and 12000x SEM Micrographs of 746 nm Wire 

Damage 

E. THIRD SAMPLE BATCH 

1. Batch Fabrication Details 

The third sample batch was fabricated with a focus on 
avoiding both overprocessing and direct exposure to the FIB 
etching, the presumed causes of the decreased luminescence 
in previous batches. A heterostructure sample was cut into 
five parts and in order to further test damage theories a 
pressed In protective cover was placed on one of the 
samples. All of the samples were then put in the FIB 
chamber. 

A 50 ym x 100 ym area was FIB etched at 500 pA for 0.2 
seconds and another 23 ym x 33 ym area was FIB etched at 
500 pA for 2 minutes on the In protected sample. Two ~0.5 
ym wire structures were also etched on this sample using 

both an 8 minute 500 pA FIB exposure and 40 minute 30 pA 

FIB exposure. 

Wire structures were then FIB etched into three of the 
remaining four samples, with the fourth serving as a 

control sample. The samples were then rinsed in HCl 

solution for 30 seconds and pressed In contacts were added 















to individual samples as required for application of 
external bias. 

2 . Batch Sample Dimensions 

The third sample batch was fabricated with one larger 
(~6.5 ym) and three smaller (~0.5 ym) wire structures. A 
significant increase in the consistency of the wire widths 
was also observed over the last two batches, and one 300 nm 
wire was constructed. Average sample dimensions were 
acquired through line scans from the SEM and the degree of 
uncertainty in the measurements correlates directly to the 
intensity contrast, however, in the case of the 300 nm 
wire, the resolution limitations of the SEM rather than the 
wire quality determined the errors in measurement. 



Table 4. Third 
Figure 22. 1OOx 


Sample 

Number 

Wire Width 
(nm) 

Error 

(nm) 

Error 

(%) 

1 

6529 

±35 

0.54 

2 

489 

±15 

3.07 

3 

571 

±20 

3.50 

4 

300 

±30 

10.0 


Batch Sample Dimensions 
SEM Micrograph of Third Batch Samples 
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F. FOURTH SAMPLE BATCH 

1. Batch Fabrication 

The fourth batch of wire structures was etched with a 
FIB directly into the already contacted and functional 
sample that was used as a control in the third batch of 
experiments. This was done to avoid contacting difficulties 
present in previous samples and to have extensive pre-FIB 
exposure drift and intensity measurements for comparison 
with the fourth batch samples. 

The use of box templates in the FIB patterning 
software allowed for a more accurate FIB pattern on the 
sample, eliminating the width variations occasionally 
produced by line template patterning used previously. The 
sample was etched with 10 minute, 500 pA series exposures 
of the FIB. 

2. Batch Sample Dimensions 

A series of four wires were fabricated on the fourth 
sample batch, ranging in width from approximately 0.8 to 9 
ym. Figure 23 shows a high definition lOOOx SEM micrograph 
of the fourth batch samples etched into the GaAs double 
heterostructure. The vertical dimension of the wires is 
~153.5 ym. 
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Figure 23. lOOOx SEM Micrograph of Fourth Batch Sample 


Average sample dimensions, shown in Table 5, were 
acquired through line scans from the SEM and the degree of 
uncertainty in the measurements correlates directly to 
resolution limitations of the SEM. 


Sample 

Number 

Wire Width 
(nm) 

Error 

(nm) 

Error 

(%) 

1 

9258 

±20 

0.36 

2 

4850 

±25 

0.36 

3 

1632 

±20 

1.16 

4 

792 

±20 

4.69 


Table 5. Fourth Batch Sample Dimensions 
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IV. EXPERIMENTAL RESULTS 


A. INTRODUCTION 

The results of the characterization work can be broken 
down into two main categories, intensity and drift 
measurements on the top-down GaAs wire structures and the 
initial observation of luminescence from the bottom-up GaN 
and ZnO nanowire structures. The top-down wire structure 
results are further divided into their respective batches, 
while the bottom-up structures are divided into SEM and CL 
imaging, CL spectra and CCD imaging. 

B. TOP-DOWN WIRE STRUCTURES 

1. First and Second Sample Batches 

The first and second batches represented initial 
attempts at production of samples displaying ID 
confinement. As discussed in the previous chapter, 
significant damage was observed after FIB etching and 
recorded via luminescence decreases in both the first and 
second batch samples. This was initially attributed to the 
multi-step processing that the samples underwent. 

Spot intensity measurements were done on the second 
batch samples and an overall pattern of decreasing average 
intensity was observed as the wire structure dimensions 
decreased. Average spot intensities are shown for the four 
different second batch sample wire structures under the 
same electron beam excitation of 1x10"^° A at 30 keV in 
Figure 24. 
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Log(Wire Width) (Microns) 

Figure 24. Comparison of Average Second Batch Wire Structure 

Intensity 


2. Third Sample Batch 

The third sample batch provided a quantitative 
measurement of the FIB damage to the samples, a comparison 
of wires etched at different FIB currents and initial 
contacting and drift measurements. The third batch control 
sample was successfully contacted and drift measurements 
were carried out on it, but those results are presented in 
the next section due to their relevance to batch four. 


The third batch samples were contacted using pressed 
In contacts and the resistances of the contacted samples 
were measured using a fine tipped probe station. Resistance 
values were far below the control sample scaled values, as 
displayed in Table 6. The significantly lower than 
predicted resistances and the inability to observe drift in 
the third batch samples led to concern for the quality of 
electrical contacts on the samples and raised the 
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possibility of parallel conduction paths due to surface 
damage. 


Sample 

Control 

2 4 ym 

14 ym 

1.4 ym 

0.7 ym 

Measured Resistance 

(kQ) 

0.3 

7 

9.4 

41.6 

120 

Predicted Resistance 

(kQ) 

N/A 

18.5 

31 

251 

592 


Table 6. Third Batch Contacted Sample Resistance 


Using an In protective layer on one of the samples 
allowed for three regions of different FIB exposure to be 
tested. A 50 ym x 100 ym area was FIB etched at 500 pA for 
0.2 seconds and a second 23 ym x 33 ym area was FIB etched 
at 500 pA for 2 minutes. The CCD image in Figure 25 shows 
the intensity contrast between the 0.2 and 120 second FIB 
exposed regions (light and dark respectively). The average 
intensity results on the In protected and FIB exposed 
regions of the sample are displayed in Table 7. 



Figure 25. 2000x, 30 keV, 1x10 A, 5 Second SEM Picture Mode 

CCD Exposure 
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Region of Interest 

Average Intensity (CCD Counts) 

In Protected Region 

7201 

0.2 Sec, 500 pA Exposure 

2391 

120 Sec, 500 pA Exposure 

1994 


Table 7. Average Maximum 0.25 Second Spot Intensities, 30 keV, 1x10 A 

Third batch samples also served to examine the 
intensity differences of wires etched at different FIB 
currents. Table 8 displays the average intensity from the 
489 nm 8 minute 500 pA FIB exposure and the 571 nm 40 
minute 30 pA FIB exposure wire structures. This was 
evidence that the exposure amperage in the FIB did not 
change the luminescence of the wire structures 
significantly. 


Region of Interest 

Average Intensity (CCD Counts) 

489 nm, 500 pA Exposure 

544 

571 nm, 30 pA Exposure 

584 


Table 8. Average Maximum 1 Second Spot Intensities, 30 keV, 1x10 A 

3. Fourth Sample Batch 

In order to control as many factors as possible and 
potentially observe carrier drift in the samples, the 
fourth batch of wires was etched with a FIB directly into 
the already contacted and functional sample that was used 
as a control in the third batch of experiments. Fourth 
batch sample wire structures displayed higher luminous 
intensities and greatly decreased variations in width due 
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to the refinement of the FIB etching process employed to 
fabricate them. 

Extensive spot intensity measurements were done on the 
fourth batch samples. Consistent with the second batch 
observations, as the wire structure dimensions decreased, 
an overall pattern of decreasing average intensity was 
observed. The average intensity observed on the widest 
structure is a factor of ~65 less than the initial (pre-FIB 
processed) material. The average spot intensities in the 
four fourth batch sample wire structures are shown in 
Figure 26. 


400 



Log(VVire Width) (Microns) 


Figure 26. 


Comparison of Average Fourth Batch Wire Structure 

Intensity 


In order to analyze drift in the fourth batch wire 

structures, it is important to first examine drift data 

from the contacted heterostructure used as the control 

sample in the third batch experiments. Figure 27 shows data 
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obtained from the sample before it was exposed to the FIB 
at fields of zero and ~0.01 V/ym and Figure 28 shows line 
scans of the same unexposed sample. 



Figure 27. 30 keV, 1x10 A, 0.35 Second Spot on GaAs 

Heterostructure, 0 and 20V bias (170 pm x 160 pm) 
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Figure 28. 30 keV, 1x10 A, 0.35 Second Spot Linescans on GaAs 

Heterostructure, 0 and ±0.01 V/pm bias 


Using the same contacts, potential differences of up 
to 20 volts were applied to the sample after the FIB 
etching of wire structures. A sample resistance of 908 Q 
was measured and there was optical evidence of an applied 
field. Changes in spot intensity from 152 counts per second 
at no bias to 33 counts per second at 20V bias also serve 
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as a clear indication that field was applied to the sample. 
Assuming that the voltage dropped at the contacts remained 
similar, an approximate potential difference of ~0.13 V/ym 
was applied during experiments. Figure 29 shows two CCD 
images taken at different bias on the smallest (~0.8 ym) 
wire structure and Figure 30 shows line scans across the 
same data. These results are further discussed and analyzed 
in Chapter V. 



Figure 29. 30 keV, 1x10 A, 20 Second Spot on 0.8 ym wire 

structure, 0 and 20V bias (136 ym x 117 ym) 



wire structure, 0 and ±0.13 V/ym bias 
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C. BOTTOM-UP WIRE STRUCTURES 

1. SEM and CL Micrographs 

The first step in expanding the transport imaging 

technique to nanowires involved locating the nanowires on 
their respective substrates and recording SEM and CL 
micrographs to determine the spectra and magnitude of 

luminous emission. It is important to note that the CL 

system collects all light emitted from any point in the 

image and maps it to the point on the image that is being 
exposed to the electron beam, so it is used to provide a 
general idea of the panchromatic intensity of the sample or 
the intensity of the sample at a selected wavelength, but 
its spatial information assumes that charge generation and 
luminescence arise from the same point. 

ZnO and GaN nanowires were deposited on Gold and Si02 
substrates but could only be located on the Si 02 substrate, 
possibly due to chemical interaction and decomposition of 
the nanowires with the Au substrate. ZnO and GaN nanowires 
and nanowire groups were observed to be bright on SEM 
micrographs mainly due to charging because of their 
electrical isolation on the Si02 substrate, figures 31 and 
32 show SEM micrographs of ZnO and GaN nanowire groups at 
2000x and isolated nanowires at a magnification of 4000x on 
the Si02 substrate. The measured diameters of 10-200 nm and 
aspect ratios of 10-100 of the nanowires were consistent 
with expected values [17,18,19] for both GaN and ZnO. 
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Figure 31. 2000x and 4000x SEM Micrographs of ZnO Nanowires, 15 

keV, 1x10"^° A 



Figure 32. 2000x and 4000x SEM Micrographs of GaN Nanowires, 15 

keV, 1x10"^° A 

After the initial SEM micrographs were taken, the CL 
system was activated in parallel with the SEM, allowing for 
simultaneous SEM and CL imaging. The CL system was set to 
panchromatic mode in order to achieve maximum intensity for 
the micrographs. The results for the ZnO nanowires are 
displayed in Eigure 33. 
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Figure 33. 2000x SEM (left) and CL (right) Micrographs of ZnO 

Nanowires, 15 keV, 1x10“^° A 

CL images were also taken at 2000x on the GaN 
nanowires in panchromatic mode, as shown in Figure 34. The 
panchromatic CL image was nearly featureless in the GaN 
case at magnifications of 2000x due to strong substrate 
emissions in the background; however some individually 
luminescent wires were visible. 
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Figure 34. 2000x SEM (left) and CL (right) Micrographs of GaN 

Nanowires, 15 keV, 1x10”^° A 

Further examination of the strongly luminescent 
regions on the GaN nanowire sample at higher magnification 
led to successful CL imaging on that sample. Individual GaN 
nanowires are visible in the lOOOOx CL image in Figure 35. 
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Figure 35. lOOOOx SEM (left) and CL (right) Micrographs of GaN 

Nanowires, 15 keV, 1x10“^° A 

2. Cathodolviminescence Spectra 

CL spectra were taken on the ZnO and GaN nanowire and 
substrate regions in order to determine the wavelength 
distribution of the panchromatic data taken in the CL 
micrographs. CL spectra were taken between 300 nm and 700 
nm with steps of 1 nm, to identify characteristic nanowire 
emission and distinguish it from the substrate emission. 

ZnO nanowires displayed a very strong peak at about 
378 nm (3.29 eV), visible in Figure 36. 



Figure 36. CL Spectra of ZnO Nanowires and Si02 substrate, 15 

keV, 1x10"^° A 
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The GaN nanowire spectra had to be directly compared 
to substrate spectra in order to clearly indentify an 
emission peak. In Figure 37, a small deviation from the 
substrate spectra can be seen in the GaN nanowire spectra 
in the vicinity of 350-450 nm and corresponds to a peak at 
about 380 nm (3.27 eV) . 




Figure 37. CL Spectra of GaN Nanowires and Si02 substrate, 15 

keV, 1x10'^° A 

The observed spectral peaks associated with the 
nanowires are in agreement with expected values of ZnO and 
GaN bandgaps of ~3.3 eV and ~3.4 eV respectively at 300K 
[16,17] . 

3. CCD Imaging of Nanowires 

The data gathered during CL clearly indicate that the 
nanowires luminesce; however, the CL system employs a high 
voltage photomultiplier while the CCD system records 
directly emitted intensity without any optical 
amplification. In order for the transport imaging technique 
to be employed on nanowires they must be visible on the CCD 
system. 
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The final results taken on the nanowire samples were a 
series of unfiltered CCD images in which GaN and ZnO 
nanowires' spatial luminescence was resolved with the OM 
and CCD system used for transport imaging. In Figure 38 a 
group of GaN nanowires are imaged with the CCD, in Figure 
39 two single ZnO nanowires are imaged. 



Figure 38. 2500x, 15 keV, 1x10 A, 120 Second Image of GaN 

Nanowires 
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Figure 39. 


2000x, 30 keV, 1x10 A, 120 Second Image of ZnO 
Nanowires 
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V. 


CONCLUSIONS AND SUGGESTIONS FOR FURTHER 


RESEARCH 

A. SUMMARY AND CONCLUSIONS 

1. Top-Down Wire Structures 

Throughout the research on top-down wire structures 
four separate sample batches were fabricated and various 
intensity and drift measurements were made on the samples. 
Extensive experience with the FIB system and samples 
allowed for a progressive refinement of the FIB etching 
process employed to fabricate the structures, however, 
there was a significant decrease in average luminous 
intensity in all samples that were exposed to the FIB. Many 
wire structures displayed evidence of electrical contacting 
problems and even those which were contacted effectively 
did not display measurable drift. 

Initial FIB etching resulted in average decreases in 
luminous intensity of over a factor of 1000 from the 
unexposed AlGaAs/GaAs/AlGaAs double heterostructure and 
wire structures of poorly defined edges and widths. 
Reduction in the processing of the sample and experience 
with the FIB system and the patterning software yielded 
significant increases in the wire structures' quality and 
reduced the average decreases in luminous intensity to a 
factor of 65. 

Although the latter sample batches' loss of intensity 
and structural variation was significantly lower than 
earlier batches, all the top-down wire structure batches 
demonstrated a consistent pattern: Figure 40 shows that in 
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all measured samples, as the wire structure dimensions 
decreased so did the average intensity. 
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Figure 40. Comparison of Average Second and Fourth Batch Wire 

Structure Intensities 



Changes in electrical properties of samples due to FIB 
exposure have been recorded in the past and characterized 
as ion-induced surface defects and passivation of dopants 
due to surface amorphization [20,21] . Ion induced surface 
defects could serve as extra scattering centers in the 
vicinity of the FIB etching which could change the non- 
radiative lifetime of the carriers within the material and 
subsequently, to first order, the overall lifetime and 
luminous intensity observed. This can be modeled as an 
effective decrease in diffusion length, following: 

, 1 1 1 

(18) - = —+- (19) 

\ e T 

where p is the mobility, i is the lifetime, ir is the 
radiative lifetime and inr is the non-radiative lifetime 
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associated with an increase in scattering centers and point 
defects in the material. 

Examination of the fourth top-down wire structure 
batch allowed a direct comparison between previously 
undamaged material which displayed drift and the same 
material after the FIB process. The unexposed sample 
displayed, on average ~26,000 counts per second on the CCD, 
while the same regions after indirect exposure to the FIB 
had intensities averaging ~400 counts per second, a factor 
of 65 decrease in intensity. 

Employing Equations (18) and (19) a new Ld of 0.44 ym 
can be calculated under the assumption that to first order, 
changes in intensity can be interpreted as changes in 
lifetime. Figure 41 shows a drift model using these 
modified parameters compared with the experimental results 
for smallest (~0.8 ym) fourth batch wire structure. 



Figure 41. 0 and 0.13 V/ym bias Ltj=0.44 ym, n=30 Simulation 

(left) and 30 keV, lxl0~^° A, 20 Second Spot Linescan 


While there was a clear broadening of the normalized 
spot under applied bias in both directions, there was no 
reliable evidence of directional drift as illustrated in 
the modeling. The main limitation to drift measurements on 
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this scale is the lack of optical resolution on the CCD 
system, as even modeled data displayed directional drift 
>1.5 ym, less than 4 pixels on the CCD image. 

The change in effective drift length due to a change 
in lifetime of the minority carriers produces a model which 
predicts drift below the resolution of our optical system, 
which is consistent with the observed data. This hypothesis 
also effectively predicts the observed trend of decreasing 
intensity with decreasing wire dimensions because if a 
constant area around the FIB etched areas is damaged then 
the smaller wire structures would display a decreased 
luminous intensity as a larger fraction of their area would 
have been damaged. 

Although the methods of FIB etching sample fabrication 
improved dramatically throughout the course of research and 
the loss of average luminous intensity in the samples was 
significantly reduced, we can conclude that, on the current 
apparatus, until further work has been done to characterize 
FIB damage, FIB etched top-down fabrication on GaAs double 
heterostructure material is not a viable method of imaging 
ID confinement and employing the transport imaging 
technique. Throughout the sample testing, reduction in the 
average luminous intensity and the correlated reduction in 
diffusion length yielded changes below the optical 
resolution thresholds for our experimental apparatus. 


2. Bottom-Up Wire Structures 


ZnO 

and 

GaN nanowires 

on Si 02 substrates 

were 

identified 

and 

imaged 

using the 

SEM, panchromatic CL 

and 

the CCD 

camera. 

Nanowire 

luminescence was 

also 
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characterized by taking CL spectra of the nanowires and the 
Si02 substrate. ZnO and GaN nanowires on Au substrates were 
not observable. 

Both the GaN and the ZnO wires showed spectral maxima 
that occurred in the vicinity of 380 nm, at a higher energy 
than most of the broadband spectra of the Si 02 substrate. 
Short wavepass filters could easily be employed to reduce 
the amount of background light in the system during 
measurements and further decrease the substrate intensity 
(noise) in CCD pictures in future experiments. 

Further research into the Au substrates' chemical 
interaction with the GaN and ZnO nanowires to determine the 
causes of nanowire disintegration or lack of luminescence 
could also be pursued, however this is less critical 
because future wire structures will be on Si 02 substrates. 
We can conclude that luminescence of the quantum wires can 
be imaged and that their emission spectra are consistent 
with expected values. This data serves as a strong 
indication of the feasibility of expanding the transport 
imaging process to ID bottom-up nanowires. 

B. SUGGESTIONS FOR FURTHER RESEARCH 

1. Advances in Bottom-Up Nanowire Structures 

The fabrication of contacted nanowires is the next 
step in the application of the transport imaging method to 
ID nanowires. Currently nanowires are grown on the 
substrate surface by CVD and remain uncontacted. Research 
groups fabricating single nanowire devices have reported 
the ability to move individual nanowires and contact them 
[4,18,19] . 
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Contacted nanowires would allow for initial 
qualitative observation of drift in ID structures, and 
coupled with increased system optical resolution could 
yield the first quantitative drift measurements. High 
resolution drift measurements could then be used to extract 
transport parameters directly from CCD images of the 
nanowires. 


2. Increased System Resolution for Transport Imaging 

In this work, the use of an optical microscope, 
internal to an SEM, to directly image the spatial 
recombination of charge generated by electron irradiation 
in nanowires has been demonstrated. However, the current 
system resolution is ~400 nm/pixel and even in these 
preliminary examinations of nanowire structures, the 
apparatus has been pushed near the limits of its optical 
resolution. 

A new National Science Foundation grant has been 
obtained to advance the technique to higher resolution 
using near field optical techniques. A near field 
microscope, with special open architecture for use inside a 
SEM, could increase resolutions of ~50-100 nm, and has been 
ordered from Nanonics Inc. and should be installed by 
summer 2006. 

Increased system resolution would be a very important 
step towards quantitative transport imaging measurements in 
nanowire structures as well as in materials which display 
very short diffusion lengths such as GaN. System 
resolutions of 50-100 nm would also allow for empirical 
verification of the reduced diffusion length hypothesis in 


60 



FIB etched top-down structures and for imaging of ID 
confinement in bottom-up nanowire structures. 

3. Transport Imaging as a Measurement of FIB Damage 

Most analyses of FIB damage use MonteCarlo simulations 
models and analysis of damage layer thickness and damage 
processes by TEM [20]. TEM studies have also shown that the 
damage layers are amorphized by the FIB, but results 
gathered from the top-down wire structures indicate that 
FIB etching may produce electrical damage that is more 
extensive than structural damage visible by TEM. 

An in depth exploration of EIB damage by means of 
transport imaging could provide a very sensitive 
measurement of sample damage at an electronic or point 
defect scale. Quantitative studies would have to be 
undertaken in different materials and different EIB 
exposure times and energies. 

More research and experimentation would also have to 
be done to determine if the observed decreases in luminous 
intensity do in fact correspond to decreases in lifetime 
within the material. Confirmation of the hypotheses 
presented in this paper could involve independent Time 
Resolved Photoluminescence (TRPL), in which the sample is 
excited by a picosecond pulsed laser and radiative carrier 
lifetime can be accurately measured. 
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APPENDIX A: DATA EXTRACTION PROGRAM 


1 %Data.m Made by Stephen winchell 

2 %Last Modified, Feb 3, 2006 

3 %Function to extract the Row and Col Max from CCD Image 

4 %Maxima marked by cross for visual error-check 

5 %Median noise filtering and noise subtraction from max value 

6 

7 function [colmax, rowmax]=data(pic); 

8 

9 %input image file to matlab and convert it to double to work with 
later 

10 I = imread(pic); 

11 I = double(I); 

12 imagesc(I) 

13 

14 %do a median noise filtering on the image 

15 J = medfilt2(I) ; 

16 imagesc(J) 

17 

18 %take off the noise floor by subtracting the average of the 
smallest sum of a row or col 

19 

20 [rows, cols] = size(J); 

21 colsum = sum(J); 

22 rowsum = sum(J'); 

23 

24 if (min(colsum)/cols) <= (min(rowsum)/rows) 

25 noisefloor = min(colsum)/cols; 

26 else 

27 noisefloor = min(rowsum)/rows; 

2 8 end 

29 

30 %find the maxima and plot it 

31 maxcols = max(J); 

32 maxrows = max(J'); 

[colmax, xposition] = max(maxcols); 
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33 



34 

[rowmax. 

yposition] = max(maxrows 

35 

colmax = 

colmax - noisefloor 

36 

rowmax = 

rowmax - noisefloor 

37 



38 

hold; 


39 

plot(xpo 

sition, yposition, 'c+' ) 
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